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ABSTRACT 

The presence of multiple populations in globular clusters has been well established thanks to 
high-resolution spectroscopy. It is widely accepted that distinct populations are a consequence 
of different stellar generations; intra-cluster pollution episodes are required to produce the pe¬ 
culiar chemistry observed in almost all clusters. Unfortunately, the progenitors responsible 
have left an ambiguous signature and their nature remains unresolved. To constrain the candi¬ 
date polluters, we have measured lithium and aluminium abundances in more than 180 giants 
across three systems: NGC 1904, NGC 2808, and NGC 362. The present investigation along 
with our previous analysis of M12 and M5 affords us the largest database of simultaneous 
determinations of Li and A1 abundances. Our results indicate that Li production has occurred 
in each of the three clusters. In NGC 362 we detected an M12-like behaviour, with first and 
second-generation stars sharing very similar Li abundances favouring a progenitor that is able 
to produce Li, such as AGB stars. Multiple progenitor types are possible in NGC 1904 and 
NGC 2808, as they possess both an intermediate population comparable in lithium to the first 
generation stars and also an extreme population, that is enriched in A1 but depleted in Li. A 
simple dilution model fails in reproducing this complex pattern. Finally, the internal Li varia¬ 
tion seems to suggest that the production efficiency of this element is a function of the cluster’s 
mass and metallicity - low-mass or relatively metal-rich clusters are more adept at producing 
Li. 

Key words: globular clusters: individual (NGC 1904, NGC 2808, NGC 362) -stars: abun¬ 
dances -stars: Population II 


1 INTRODUCTION 


Globular clusters (GCs) display internal variations in elements af- 
fect ed by proton-capture reactions (e.g., C, N, O, F, Na, Mg, and 
Al; iGratton, Carretta & Bragaglial 2012ti and exhibit spl it evolu¬ 
tionary sequences l lPiotto et al.ll2012l : iMilone et alj[20l^ in their 
colour-magnitude diagrams. The presence of both features requires 
that GCs comprise multiple populations and it is this revelation 
that has renewed interest in these old stellar systems. In order to 
decipher their complex chemical history, a concerted and system¬ 
atic effort has been made to survey GCs for key diagnostic chem- 


* Based on observations taken with ESO telescopes under program 094.D- 
0363(A) 


ical species. Abundances have been determined for a considerable 
number of stars, from the main sequence through to the asymp¬ 
totic g iant branches (AGB) , via dif f erent approaches and techniques 


(e.g., 

Carretta et al. 200fil. 20071; Johnson & Pilachowski 20 id; 

Gratton et al. 

2011 

I 2 OI. 1 2014;ICamDbell et alJ2013 

: Marino et al 

2011, 

2013, 

2014 

; Yone et al, 20131. 

2014 

2 OI. 5 I: 

Johnson et al 

2015; 

Lamb et alj 

2 OI 5 I: Valenti et al.l 

2015 

Cordero et al, 20151; 

Mucciarelli et alJl2015L and references therein). 


Red giant branches have been extensively studied with the 
largest survey of chemical compos itions for GC giants under¬ 
taken by Carretta and collaborators dCarretta et alj|200^ l2009aL 
I 2 OI lLl20r3l . l2014h . They determined iron-peak, or-, proton-capture, 
and neutron-capture element abundances for thousands of giants 
in 25 GCs, characterised by different morphological properties 
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(i.e., mass, metallicity, age, environment, concentration, horizontal 
branch morphology). Recently, the northern sky has been surveyed 
by the SDSS-111 Apache Point Observatory Gala ctic Evolution Ex¬ 
periment (APOGEE) and iMeszaros et all JlOlSh provided light el¬ 
ement abundances for approximately 400 giants across ten GCs. 

To date, almost every system (see how ever the peculiar cases 
of Ruprecht 10 6, Terzan 7, and Pal 12; IVillanova et alj l2013l ; 
ISbordone et alJl200^ ; ICohenll2004h analysed displays evidence for 
multiple stellar populations, suggesting that this feature is in fact 
ubiquitous. Spectroscopically this evidence is in the form of the 
characteristic GC abundance patterns which are present in the con¬ 
stituent dwarfs and giants and therefore not specific to environ¬ 
ment or evolutionary stage. The observed abundance patterns are 
expected in material that has undergone proton captures at high 
temperatures - enhancements in N, Na, and Al are accompanied 
by corresponding depletions in C, O, and Mg abundance^. The 
C-N and Na-0 anticorrelations display very different levels of en¬ 
hancement/depletion across different systems suggesting that the 
pollution mechanisms and the nature of the polluters may actually 
vary from cluster to cluster. Observations of the Mg-Al anticorre¬ 
lation (although not always present) imply a similar conclusion; 
the multiple populatio ns are universal, but the progenitors are not. 
ICarretta et al.l ( l2009al) found that the extent of the Al enrichment 
(and possible corresponding Mg depletion), is not the same in all 
GCs - it i s related to a combin ation of cluster mass and metal- 
licity. The Carretta et alj ( l2009all result has been corroborated by 
iMeszaros et alj 1 2015h for a different set of systems. They con¬ 
firmed that the spread in Al abundances increases significantly as 
cluster average metallicity decreases. They interpreted this as evi¬ 
dence that low-metallicity, intermediate-mass AGB polluters were 
more common in the metal-poor clusters. In add ition to AGB 
stars l lCottrell & Da Co^ll 98 ll ; I Ventura et alj|200lh . other candi¬ 
dates suggested as viable sources of i nternal pollution includ e fast- 
rotating mass ive stars (FRMS, e. g., IPecressin et alj|2 007h. mas¬ 
sive b inaries jde Mink et alj[200^ . or novae ( Maccarone & Zur^ 
[2 oT 3) . Unfortunately, no one theory can satisfactorily account for 
the observed chemical patterns. 

We have commenced a large observational survey with the sci¬ 
entific motivation to shed light on the nature of the internal polluters 
in GCs. Our focus has been on the element lithium which we con¬ 
sider to be one of the most powerful diagnostics of the progenitors 
responsible. The s cientific context has been thoroughly discussed 
in our recent work jD’Orazi et alj201^ hereinafter Paper i), where 
we presented Li and Al abundances in GC giants of NGC 6218 
(M 12) and NGC 5904 (M 5). In essence. Paper i is based on a 
straightforward nucleosynthetic prediction: if the polluters produce 
no Li, then an anticorrelation should develop between Li and Al. 
This is because the second generation stars (Al-rich) will form from 
(almost) Li-free ejecta, whilst the primordial population are char¬ 
acterised by high Li abundances (appropriate to a value initially 
around the value of the Spite plateau, but modified by first dredge- 
up -FDU) and a modest Al content (as field stars at the same cluster 
metallicity). Furthermore, this anticorrelation would be erased if Li 
production occurs across the different stellar generations. Evidence 
for Li production in the progenitors was previously detected in the 


GC M4 bv IP’Orazi & Marini 120 id) and sub sequently confirmed 
by other studies (e.g., Mucciarelli et al.ir201 ih . In the two clusters 
analysed in Paper i, we determined that Li production in the progen¬ 
itors is necessary to explain the class of stars that are both Li and Al 
rich. The conclusion favours AGB stars as responsible for internal 
pollution because their chemical yie lds ca n predict Li production 
via the so called Cameron & Fowler ( 1 197 ih mechanism. Curiously, 
the Li content measured in second-generation stars is exactly the 
same as the primordial population, requiring a kind of ad hoc fine 
tuning in the Li synthesis process. Current stellar theory does not 
allow for such Li production in other candidates. 

All three clusters. Ml2, M4 and M5 are of similar metallicity, 
with the last the most massive. In M5 we identified an abundance 
pattern not present in M4 and Ml2: the cluster hosts an extreme 
population of stars that are Li-deficient and Al-rich but which can¬ 
not be explained by a simple dilution model. Moreover, our com¬ 
parison with previous studies yielded a correlation between the in¬ 
ternal Li spread and the current cluster mass: the Li production 
seems less efficient in more massive GCs, whereas in smaller sys¬ 
tems first-generation (EG) and second-generation (SG) stars share 
an almost identical Li content. We noted that small statistics and 
possible systematics may impact upon the speculative conclusions 
drawn from the comparison and thus we improve upon the sample 
size here. 

In the present paper, we report the results from our study of Li 
and Al abundances for cluster giants in three more clusters, namely 
NGC 1904, NGC 2808, and NGC 362. This study complements 
our previous work focussed on M12 and M5 and offers the largest 
Li survey in GC RGB stars available to date. 


2 OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

Observations were c arried out in visitor mode with FLAMES 
jPasQuini et al.|[2003) mounted at UT2 of ESO-VLT on December 
11-13 2014 (Program 094.D-0363(A)). Spectra were collected with 
the scientific aim of inferring Li and Al abundances for large sam¬ 
ples of RGB stars, below and above the RGB bump luminosity, 
in several GCs, as part of our dedicated survey (see Paper i). We 
utilised the HR15N high-resolution setup, covering simultaneously 
H„ and the Li doublet at 6708 A, which grants a nominal resolution 
of R= 17,000 and a spectral coverage from 6470 to 6790 A. We se- 
le cted our target stars ba sed on photometric information published 
in iMomanv et al.l ( l2004tl . that is stars lying on the ridge line of the 
giant branch and with no companion brighter than 2 mag within 
2 . The colour-magnitude diagrams (BV photometry) for the three 
target clusters ar e shown in Fig. [J where the location of the bump 
is marked (from l^taf et alj[2013l) . The initial samples comprised 
55, 106 and 78 giants in NGC 1904, NGC 2808, and NGC 362, 
respectively. 

Data were reduced using the dedicated 
ESO pipeline (version 2.12.2, available at 
http: //WWW. eso. org/sci/software/pipelines/), provid¬ 
ing bias-subtracted, flat-fielded, optimal extracted and wavelength 
calibrated one dimensional spectra. Continuum normalisation, 
shifting to rest frame, and combination of multiple exposures 
were then performed using IRAlj| tasks. After radial velocity 


* In the most widely accepted scenario, the fact that GCs display intrin¬ 
sic variations in their light-element content points to different stellar' gen¬ 
erations, whereby a fraction of first generation stars polluted the interstel¬ 
lar matter from which secon d generations have b orn. However, a different 
model has been proposed bv iBastian et al.l (201^. 


^ IRAF is the Image Reduction and Analysis Facility, a general purpose 
software system for the reduction and analysis of astronomical data. IRAF 
is written and supported by National Optical Astronomy Observatories 
(NOAO) in Tucson, Arizona. 
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Figure 1. Colour-magnitude diagram for the three globular clusters. The location of the RGB bump is marked ffrom l^taf et alj2013h . Stars analysed in the 
present study are emboldened. 





computations, six stars in NGC 1904, 31 stars in NGC 2808, 
and six stars in NGC 362 turned out to be not members of 
their respective clusters and were thus discarded. We obtained 
average values of ¥„<,= 205.78±0.54 km s^’ (NGC 1904), 
V„<,= 102.79±0.97 km s-‘ (NGC 2808), and 222.95±0.61 
km s“* (NGC 362 ), whi ch agree very well with measurements 
reported in Harris ill 9961 . 2010 update), that is Vrad^ 205.8±0.4 
km s“*, yrad- 101.6±0.7 km s“*, and yrad- 223.5±0.5 km s“*, 
respectively. Our results ar e also in very good agreeme n t with 
previous measurements by ICarretta et al.l l2009bl . l2006l . l2013l) 
of V™j= 205.60±0.48 for NGC 1904, V,„rf=102.90±1.01 for 
NGC 2808 and V™rf=223.31±0.53 for NGC 362. The heliocentric 
radial velocity distributions for cluster members are shown in 
Figure [2 


As for Paper i, stellar parameters were derived in the fol¬ 
lowing manner. Initial effective temperatures (T^g) were calculated 
from {V - Ks) colo urs (with optical photometry retrieved from 
iMomanv et al.l2003 and infra-red magnitud es from the 2MASS 
poi nt source catalogue, SkrutsWe^et^ 2006) and the calibration 
by lAlonso. Arribas & Martmez-RogeJ 1 1999 ). To do this we as¬ 
sumed reddening and metallicity values collected in the Harris’ 
catalogue, that is E(B-V)=0.01 and [Fe/H]=-1.60 for NGC 1904, 
E(B-V)=0.22 and [Ee/H]=-1.14 for NGC 2808, and E(B-V)=0.05 
and [Ee/H]=-1.26 for NGC 362. As done in our previous stud¬ 
ies, we co nverted E(B-V) values to E(V-K) adopting the rela¬ 
tionship bv ICardelli. Clayton & Mathis ( Il989h . The final, adopted 
temperatures come from a relationship between those T^g and V 
magnitudes. Our sam ples consist of 22 RGBs in common with 
ICa^etta^ ai] (l2009bt) for NGC 1904 and 21 stars in NGC 362 al¬ 
ready presented by Carretta et all ( l2013l) . the mean difference in 
temperature is A7’eff=5.95±5.90 K and A7’cff=-32.48±5.00 K. The 
larger difference for NGC 362 is probably related to the fact that 
in Carretta’s study ( l2013h . T^g values come from the relationship 
between Ks magnitudes and VK col ours, whereas we use d the V 
magnitude (as done for NGC 1904 bv ICarretta et alj|2009bl) . 


Surface gravities (log g) were calculated assuming 


those Tqg, M;,„;j;;=4.75, the bolom etric correction from 
lAlonso, Arribas & Martmez-Rogeil h999l) . a mass of M=0.85Mo, 
and distance moduli from Harris, by using the following equation: 


log — = log — 

go Ma 


■log — + 41og 

Lq 


Ttg 

Feff.G 


Einally, mi croturbulence velocities were comp uted with the 
relationship by ICratton, Carretta & Castellil l ll996li . that is f = 
2.22 - 0.322 X log g 

Similarly to Paper i, we gained abundances for Li and Al from 
the strong doublets at 6708 A and 6696/6698 A, respectively. Due 
to the occurence of cosmic rays on the spectral features under con¬ 
sideration and/or to lower S/N ratios, we could obtain Li abun¬ 
dances for 47 stars in NGC 1904, 68 stars in NGC 2808, and 67 
stars in NGC 362. The Al content has been inferred for 48, 65, 
and 66 cluster members, respectively. We performed spectral syn¬ 
thesis calculations using the local t hermodynami cal equilibrium 
(LTE) code MOOG by C Sned en ( ISnedenfll97i 2014 version) 
and the lCastelli & Kurucj (|2004|) grids of model atmospher^ with 
tt-enhancement (-1-0.4 dex) and no convective overshooting. An ex¬ 
ample of spectral synthesis for the Li feature is displayed in Fig¬ 
ure [3] We finally applied NLTE corrections to our Li abundances 
following prescriptions bv iLind et alj ( l2009l) . 

Random uncertainties were calculated as described in Paper i, 
to which we refer the reader for an extensive discussion. Here we 
briefly recall that (internal) error measurements include two kind 
of contributions: uncertainties related to the best-fit determination 
and errors due to the adopted set of atmospheric parameters. Given 
their independence, they can be added in quadrature in order to 
derive the final internal uncertainty on Li and Al abundances (see 
Paper i for details). Typical values range from 0.08 - 0.12 for Li and 
0.10-0.15 for Al. 


^ Available at http: //kurucz. harvard. edu/grids. html 
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Figure 2. Radial velocity distributions for our target stars. 
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Figure 3. Examples of spectral synthesis for one star in each target cluster. 


3 RESULTS AND DISCUSSION 

Results of our analysis are reported in Table|T] optical and 2MASS 
Ks photometry, heliocentric radial velocity, S/N ratios as calculated 
around the Li i doublet and final stellar parameters (effective tem¬ 
perature, gravity and microturbulence) are listed for each target star 
in the three clusters (the full table is available in the online version). 
Both the LTE and NLTE Li abundances are presented in Columns 
10 and 11, respectively, while [Al/Fe] ratios are given in Column 
12 (asterisks denote upper limits in their respective measurements). 


3.1 The trend of lithium with luminosity 

We plot Li as a function of V magnitude in Figure|4]for NGC 1904 
(left-hand panel), NGC 2808 (middle panel), and NGC 362 (right 
panel). Th e locations of the RGB bump in each cluster are 
taken from iNataf et ^ ( 1201 3h and are marked by vertical lines 
at Vi,„,„p=15.88, 16.24 and 15.40, respectively. As expected from 
standard stellar evolution, stars that have already experienced the 
FDU are depleted in Li by factors 15-20 with respect to the Spite 
value (thus if we assume that this is the pristine Li abundance, 
with respect to their original content, although some further ex¬ 
tra mechanisms such as atomic diffusion might have played a mi¬ 
nor role). We determined mean lithium abundances of A(Li)NLTE= 
0.97 ± 0.02 (rms=0.11, 23 giants) for NGC 1904, A(Li)NLTE= 1-06 
± 0.02 (rms=0.13, 65 stars) for NGC 2808 and A(Li)NLTE=L02 


± 0.01 (rms=0.09, 57 stars) for NGC 362 by considering only 
those stars fainter than the bump luminosity. All three clusters share 
the same mean Li abundance once observational uncertainties are 
taken into account, and furthermore, they agree with the Li abun¬ 
dances determined for every other GC giant currently in the liter- 
atur^ If we apply NLTE corrections to data published in Paper i, 
we obtain A(Li)NLTE=1.07±0.01 for NGC 6218 and L02±0.01 for 
NGC 5904 (average abundances under the local thermodynamical 
equilibrium approximation were 0.98±0.01 and 0.93±0.01, respec¬ 
tively). We stress that giants for these five GCs are analysed in 
a strictly homogeneous way (same code, line-lists, model atmo¬ 
spheres and techniques) and consequently systematics do not im¬ 
pact on our conclusions. We find als o that our derived Li ab un- 
dances are similar to measurements bv lMucciarelli et aP ( 1201 4h in 
the massive GC M54. M54 is located in the centre of the S agit- 
tarius dwarf spheroidal galaxy dlbata, Gilmore & Irwin|[l99^ and 
is kn own to exhibit a sprea d in metallicity and x-process elements 
('e.g.. lCarretta et ahlbOlOall . Mucciarelli and collaborators derived 
a mean value of A(Li)=0.93± 0.11 (or 1.01, depending on the 
adopted NLTE corrections). Thus our results seem to indicate that 


'* A very small fraction, below the 1%, of Li-rich g iants are also detected in 
GCs (see e.g., our recent work. lp^Orazi et al.l2015l and references therein). 
The nature of these objects is however related to a peculiar evolution and 
hence its discussion is not relevant to the present work. 
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GC giants (before they reach the RGB bump) share the same pri¬ 
mordial Li abundances, independent of cluster’s mass, metallicity, 
horizontal branch morphology, or location in the Galaxy (disk vs 
halo) or beyond (Milky Way vi dwarf Spheroidal GCs). We note 
that it is also possible that the stars are bom with a different Li con¬ 
tent but FDU or some rapid pre main-sequence process may operate 
in such a way that they eventually end up with the same Li abun¬ 
dances. Stellar models currently do not lend support to the FDU 
scenario whilst PMS depletion may require some yet to be included 
physical process. We summarise the GC global structural parame¬ 
ters along with Li abundances in Table where we report values 
published in Harris’s catalogue for metallicity ([Fe/H]), absolute 
cluster magnitude (a proxy for the current mass). King-model cen¬ 
tral concentration (a c denotes core-c ollapsed cluster). We include 
also the horizontal branch ratios from iGratton et al.] ( I 2 OIOI 1 . the in¬ 
terquartile range of Na-O distributions as well as fractions of pri- 
mordial (P), intermed iate (I) and extreme (E) population^ from 
ICarretta et al.l ( 1201 Obh and lCarretta et alj ( 1201 3l . see also following 
discussion). 

Once stars evolve past the bump in the luminosity function, 
their Li content is efficiently depleted due to the onset of some kind 
of extra mixing processes: the H-buming shell erases the molecu¬ 
lar weight discontinuity (the so called “yu-barrier”) established by 
first dredge-up and deep mixing occurs. The nature of this kind of 
extra mixing is one of the hottest topic in current stellar evolution 
theory and different mec hanisms have been suggested such as e .g., 
“thermohaline mixing]]_(lEg gletom_^arbom & Lattanzidl2OO0 or 
magnetic buoyancy dBusso et~n .l2007h . For clusters NGC 1904 and 
NGC 362, where we well sample this region of the HR diagram, the 
location of the photometric and spectroscopic RGB bumps match 
very well: the onset of lithium depletion coincides with the bump 
magnitude. In the case of NGC 2808, their is a dearth of observa¬ 
tions at this magnitude due to an observational bias in our target 
selection. Our sampling is somewhat discontinuous but it is clear 
that stars beyond the bump luminosity possess undetectable levels 
of Li. 


3.2 Aluminium abundances 

We determined Al for 48, 65, and 66 giants (below and above the 
RGB bump) in NGC 1904, NGC 2808, and NGC 362, respec¬ 
tively (179 GC members in total). The [Al/Fe] ratios are shown 
as a function of for our target stars in Figure [5j no obvi¬ 
ous trend is evident from the data, corroborating the reliability of 
our spectral analysis. In NGC 1904, we found an average value 
of [Al/Fe]=0.19±0.05 (rms=0.35), with abundances ranging from 
-0.30 to 0.70 dex, corresponding to a variation of a factor of 10 
within this cluster. Our Al abundance determinations for NGC 1904 
are offset by approxi mately 0.45 dex to the values determined by 
ICarretta et al.l ( l2009al) . i.e., [Al/Fe]=0.64±0.15, rms=0.4L We have 
no stars in common with their sample so a direct comparison aimed 
at revealing the cause of the discrepancy in not possible. We re¬ 
call that Carretta’s work exploited high-resolution UVES spectra 
and was then focussed on brighter giants, thus NLTE effects might 
contribute to the difference. Their sample was also smaller than 
ours with only one of their eight RGB stars labelled ‘primordial’ 

^ The cluster populations ai'e defined according to their O and Na abun¬ 
dances: P stars ai'e those having [Na/Fe]<[Na/Fe]n]i„ -I- 4cr, whereas second- 
generation stars are classified as intermediate or extreme if [O/Na] is larger 
or smaller than -0.9 dex, respectively. 


([Al/Ee]=-0.216). Their remaining seven stars display different 
levels of enhancement and are consi dered second generatio n mem¬ 
bers. It is interesting to note that the lCarretta et alj ( l2009ah results 
imply that the relative distribution of different stellar populations 
determined from Al enhancement does not reproduce t he distr ibu- 
tion inferred from Na and O abundances. ICarretta et alj j2009bh ob¬ 
tained fractions of FG=40% and SG=60% (50% I + 10% E popula¬ 
tions) which better reflects our Al distributions (see Figure]^ where 
the two peaks in Al abundances have almost the sa me size). It is 
temptin g to speculate that limited sample size in the lCarretta et alj 
( l2009tJl survey may skew the statistics; however the internal vari¬ 
ation in the Al content agrees very well between the two studies 
(with Carretta’s study providing a lower limit to the Al spread of > 
0.85 dex). 

For NGC 2808 we derived a mean Al abundance of 
[Al/Fe1 =0.32±0.05 (rms=0.39), which is consistent with Carretta’s 
( l2009ah determination of [Al/Fe] =0.40±0.14 (rms=0.49) and with 
the recent analysis by ICarrettal (l2014h . i.e., [Al/Fe] =0.46±0.09 
(rms=0.47). We detected an internal spread of mo re than 1 dex (- 
0.1 to 0.98 dex) agreeing with previous findings bv iBragaglia et alj 
( l2010h based on unevolved stars. They published the first measure¬ 
ment of Al abundances for main-sequence stars in NGC 2808 using 
A-shooter@VLT spectra, deriving [Al/Fe]=-0.2 for the red main- 
sequence star (i.e., it belongs to the primordial population of the 
cluster) and [Al/Fe]=-1-1.1 for t he (blue) second -generation dwarf. 
We also confirm the finding by ICarrettal ( l2014h that the Al distri¬ 
bution is clearly trimodal in NGC 2808 (see Figure]^. A possible 
hint for a double-peak trend in the Al content seems to be present 
for NGC 1904, but the pattern is less clear (we also note that due to 
the lower resolution and signal-to-noise (S/N) ratios of our spectra, 
observational uncertainties might be responsible for smearing off 
the distribution). 

The spread of [Al/Fe] in NGC 362 spans 0.6 dex (range 
-0.25 dex to -t-0.35 dex), with an average abundance of 
[Al/Fe]=-0.02±0.02 (rms = 0.14). We find a ~ 0.25 dex offset 
compared to the value reported by bv ICarretta et alj ( l2013h . namely 
[Al/Fe]=0.24±0.05, rms=0.19. As we speculated for NGC 1904, 
the differences between the two studies might be related to sample 
size, small differences in temperature scales and/or NLTE effects 
and we are once again in the unfortunate position where our surveys 
contained no stars in common. There is yet again good agreement 
between the observed internal spreads with Carretta et al. quoting 
a peak-to-peak variation of 0.56 dex. 

Although for two of our clusters we derive different mean 
[Al/Fe] abundances to those previously published in the literature, 
our discussions and conclusions remain unaffected. For our scien¬ 
tific analysis, zero-point offsets will not impact upon the fact the 
large spread in Al is evidence for multiple populations. In Fig¬ 
ure 7 we plot o ur [Al/Fe] as a functi on o f [Na/Fe] and [0/Fe] 
as determined bv ICarretta et al.l ( l2009bh and ICarretta et al J ( l2013h 
for NGC 1904 and NGC 362 (we have no stars in common for 
NGC 2808). We still recover the pattern of a very extended positive 
correlation between Na and Al for NGC 1904 as well as a sig¬ 
nificant O-Al anticorrelation. The Pearson correlation coefficient is 
r=0.84 (18 degrees of freedom, significance level more than 99.9%) 
and r=-0.83 (16 degrees of freedom, significance level more than 
99.9%), respectively for Na a nd O. NGC 362 exh ibits a less pro¬ 
nounced Al variation (see also lCarretta et alj|2013h with respect to 
the more massive GC NGC 2808 or the more metal-deficient clus¬ 
ter NGC 1904. However, statistically meaningful (anti)correlations 
with O and Na do appear, with Pearson correlation coefficients of 
r=0.83 for Na and r=-0.71 for O. The probabilities that those corre- 
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Table 1. Photometry, stellar parameters, kinematics and chemical information for our sample stars. The complete table is made available on line. 


Star ID 

RA 

(hh:mm:ss) 

DEC 

B 

(mag) 

V 

(mag) 

K., 

(mag) 

Vrad 

(kms-‘) 

S/N 

(@6708A) 

Teff 

(K) 

log g 
(cm s“^) 

(kms-‘) 

A(Li)LTE 

dex 

A(Li)NLTE 

dex 

[Al/Fe] 

dex 

1904-4196 

05:23:58.912 

-24:30:02.012 

16.798 

16.003 

13.698 

209.82 

188 

4841 

2.20 

1.51 

0.85 

0.95 

-0.20 

1904-8925 

05:24:17.168 

-24:31:24.698 

16.858 

16.083 

13.896 

203.43 

175 

4858 

2.24 

1.50 

0.80 

0.89 

0.60 

1904-6388 

05:24:11.917 

-24:32:37.674 

16.900 

16.105 

13.916 

202.70 

180 

4863 

2.25 

1.50 

1.00 

1.09 

-0.15 

1904-7558 

05:24:13.986 

-24:32:01.426 

16.911 

16.107 

13.913 

204.79 

160 

4863 

2.25 

1.50 

0.95 

1.04 

-0.10 

1904-8251 

05:24:17.435 

-24:31:42.530 

16.912 

16.131 

13.880 

208.93 

170 

4868 

2.26 

1.49 



0.60 


hJ 
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1 1 1 1 1 1 1 1 
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Figure 4. NLTE Li abundances as a function of V magnitudes. The vertical bar marks the RGB bump. Empty circles denote measurements and upside-down 
triangles are for upper limits 


lations happened by chance are also in this case less than 0.1% . We 
thus corroborate prev ious results bv ishetrone & Keand l l2000tl and 
ICarretta et al.l ( 1201 3h who compared the second-parameter pairs 
NGC 288 and NGC 362. The light-element variations in NGC 362 
are more extreme (a factor of two in [Al/Fe]) probably owing to 
its larger mass (they have the same metallicity but NGC 362 is five 
times larger than NGC 288). 

Our findings confirm that the extent of internal variations in 
proton-capture elements, and in particular the presence of species 
requiring very high burning temperatures such as e.g., aluminium, 
depends upon a combination of the cluster’s mass an d metallicity. 
This result was first reporte d by Carrettaet aUj^OPah and later on 
verified in other clusters bv iMeszaros et alj ( l2015h . 


3.3 Lithium abundances and the multiple population 
framework 

The determination of Li abundances in conjunction with p-capture 
elements (e.g., Na, O) allows us to put strong observational con¬ 
straints on the internal pollution mechanisms and their correspond¬ 
ing stellar source. All three clusters analysed in this study exhibit 
significant internal variations in their Al content (see Section lT^ . 
implying that temperatures in excess of 70x10® K have been 
reached within the first generation stars that provided the enriched 
ejecta to give birth to subsequent stellar generations. It is quite rea¬ 
sonable to expect that at those high temperature all Li is destroyed 
(Li starts burning at only 2.5x10® K) and that their second genera¬ 
tion stars are born of Li-free ejecta (though a mix with some pristine 
material has to be postulated regardless of the polluters’ type, e.g. , 
Prantzos . Charbonnel & Iliadi dboOTl : iD’Ercole et alj[200a 1201 ll : 

Conroy & Spergelll201 ll) . It follows that the primordial GC popu¬ 


lation should possess lithium at the Spite plateau level whereas the 
second-generation stars should demonstrate lower Li abundances 
(the extreme population should have virtually no Li if we assume 
they come from pure ejecta). Our observations can test this scenario 
because our analysis would detect conventional (anti)correlations 
between Li and proton-capture elements (that is the positive corre¬ 
lations Li-0, Li-C, and Li-Mg along with anticorrelations between 
Li and Na, N, Al). If Li production took place within the polluters 
the anticorrelations are substantially erased. 

During our previous investigations, we found that Li must 
have been produced in th e progenitors that polluted the GCs M4 
dP’Orazi & Marinal2010h and Ml2 (Paper i): despite having their 
O content reduced by factors 3-4, second-generation stars in these 
clusters are still L i rich. Similar findings have been also obtained by 
iMucciarelli et alj ( 1201 4h for the extra-galactic GC M54, where Li 
abundances in the first and second generation stars are essentially 
the same. In Figure[^we present our Li abundances as a function of 
[Al/Fe] ratios (for stars below the bump luminosity). In NGC 1904 
we detected evidence for a Li-Al anticorrelation, with the Pear¬ 
son correlation coefficient r=-0.67 (significant at more than 99.9% 
level). However, the observed pattern reproduces the trend we dis¬ 
covered in the GC NGC 5904 (M5) (see Paper i): both first and 
intermediate generation stars contain the same Li abundance sug¬ 
gesting internal polluters able to produce Li, such as intermediate- 
mass AGB stars. The clusters M5 and NGC 1904 differ from M4 
and M12 in that they host an extreme population that exhibit a Li- 
deficient pattern. Since we could only derive upper limits, we are 
not able to ascertain whether this extreme population is totally Li- 
free or a small amount of Li has been produced. The other key fea¬ 
ture emerging from our study is that a simple dilution model would 
fail in reproducing the observed pattern, by assuming that the Inter- 
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Table 2. Structural properties and Li abundances (for giants below the bump luminosity) for clusters from our survey 


Cluster 

[Fe/H] 

(dex) 

M,/ 

(mag) 

c 

HBR 

IQR[Na/0] 

P 

(%) 

I 

(%) 

E 

(%) 

<A(Li)> 

(dex) 

rms 

(dex) 

NGC i904 

-1.60 

-7.86 

1.70c 

0.89 

0.759 

40 

50 

10 

0.97 

0.11 

NGC 6218 

-1.37 

-7.31 

1.34 

0.97 

0.863 

24 

73 

3 

1.07 

0.06 

NGC 5904 

-1.29 

-8.81 

1.73 

0.31 

0.741 

27 

66 

7 

1.02 

0.11 

NGC 362 

-1.26 

-8.43 

1.76c 

-0.87 

0.644 

22 

75 

3 

1.02 

0.09 

NGC 2808 

-1.14 

-9.39 

1.56 

-0.49 

0.999 

50 

32 

18 

1.06 

0.13 


CD 

pin 


<1 







Figure 5. [Al/Fe] ratios versus effective temperatures (T^g) for the whole sample. As for Figure ?? upside-down triangles denote upper limits 


mediate population is a simple mixture of primordial and extreme 
population (as found for M5, see Paper i for further details). 

The massive cluster NGC 2808 displays a very similar pat¬ 
tern in terms of Li and A1 abundances to M5 and NGC 1904: 
the correlation coefficient between Li and A1 is r=-0.57 (62 de¬ 
grees of freedom, significance level at more that 99.9%). Analo¬ 
gous to NGC 1904, the intermediate population is still significantly 
Li-rich (despite being enhanced in Al), whereas the extreme pop¬ 
ulation shows evidence of Li depletion (to different levels). It is 
noteworthy that in this case the extreme component is more numer¬ 
ous than for NGC 1904 and remains in agreement w ith the relative 
fraction s of different stellar populations derived by ICarretta et alj 
( l2009bl) according to Na and O abundances. The fraction of E stars 
in NGC 2808 is almost twice that of NGC 1904, probably because 
of its large mass. Thus, in both these GCs we require a Li produc¬ 
tion between the first and the intermediate second-generation stars 
(and hence AGB stars seem to provide us with a culprit), whereas 
we are not able to discard any of the candidate polluters to account 
for the observed pattern in extreme stars. They indeed could be ei¬ 
ther fast-rotating massive stars or asymptotic giant branch stars that 
did not efficiently produce Li (likely due to a combination of their 
mass and metallicity regimes). 

Conversely, NGC 362 displays a similar behaviour to that 
which we have observed in M12 and M4. There is no Li differ¬ 
ence across the stellar generations. A Pearson correlation coeffi¬ 
cient here of r=-0.04 is not statistically meaningful, implying no 
correlation between [Al/Fe] and A(Li) exists. In Figure [8] there are 
a handful of stars showing a slightly lower Li abundance (marked 
as asterisks in the plot). However all these giants are located very 
close to the bump luminosity, namel y they are all withi n 0.3 mag 
from the bump magnitude quoted bv iNataf et alj l l2013tl . This im¬ 


plies that their lower Li content might be actually due to evolution¬ 
ary effects and not necessarily related to the multiple populations. 
Alternatively, they could be binary stars, showing a Li-depleted na¬ 
ture as a consequence of their binary evolution. However, an offset 
between our photometry and the study by Nataf and collaborators 
is present and we indeed detect the RGB bump at 0.45 mag fainter 
than theirs (Vi,„„,p=15.46). Had we adopted this value for the bump 
luminosity, we could have discarded those stars from the present 
discussion. As a consequence, our results tentatively indicate that 
significant Li production has occurred in NGC 362 with AGB stars 
seemingly responsible. 

In Paper i we reported the detection of a negative correlation 
between the absolute cluster magnitude (a proxy for the current 
mass) and the internal Li variation (as given by the peak-to-peak 
variation, see Figure 14 in that paper). We interpreted this as evi¬ 
dence of different efficiency in the Li synthesis mechanisms within 
different clusters, with smaller GCs exhibiting a more conspicuous 
level of Li production (that makes first and second generations in¬ 
distinguishable in terms of their Li content). The relationship was 
proposed with the caveat that we were unable to discern the role 
of metallicity in such a limited sample of clusters. The results pre¬ 
sented here do in fact suggest that metallicity plays a role in de¬ 
termining the internal spread in lithium. NGC 362 is much less 
massive than NGC 2808 and has a metal content more than a fac¬ 
tor of two higher than NGC 1904 which indicates that Li produc- 
tion is less efficient in more massive and more metal-poor systems. 
ICarretta et alj ( l2009al) reached the same conclusion about the level 
of Al spread in their cluster sample. 

Finally, it would be tempting to speculate that the reason for 
the high Li content in second-generation stars might be related to 
their He-rich nature. It is possible that FDU is less efficient in these 






















8 V. D’Orazi et al. 




[Al/Fe] 



Figure 6. Histograms of Al abundances for all giants under consideration 





[O/Fe] 


Carretta + 2009 


[O/Fe] 


Carretta + 2013 


Figure 7. Run of [Al/Fe] with [Na/Fe] and [O/Fe] for stars in common with ICarretta et al] j2009tj] . Filled circles and empty triangles are for measurements 
and upper limits, respectively. No stars are in common for cluster NGC 2808. 


stars resulting in less Li depletion. We ran stellar models (Angelou 
et al., in prep) with calculations for FG and an extreme popula¬ 
tion stars which indicate differences at the level of A(Li)~ 0.1 dex 
owing to the respective convective envelope penetration. These dif¬ 
ferences lie within the observational uncertainty and are therefore 
indistinguishable. Nevertheless, if the He abundance is responsible 
we should expect high Li in E stars compared to I stars as E stars 


are thought to contain a very high level of He enrichment (up to 
Y~0.4 suggested for NGC 2808). We demonstrated that, whenever 
present, the E component is actually more Li-poor (possibly even 
Li-free) than the intermediate population. 
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[Al/Fe] 


Figure 8. NLTE Li abundances as a function of [Al/Fe] ratios for stars fainter than the bump luminosity (see text for discussion). 


4 SUMMARY AND CONCLUDING REMARKS 


The aim of our survey has been to probe the stellar sources of inter¬ 
nal pollution in GCs. We have determined abundances for Li and 
A1 in a sample of roughly 180 giant stars in clusters NGC 1904, 
NGC 2808, and NGC 362. This investigation complements our pre¬ 
vious work which focussed on NGC 6218 (Ml2) and NGC5904 
(M5) and provides the community with the largest, homogeneous 
database of GC lithium abundances available in the literature to 
date. The robustness of lithium as diagnostics of the enrichment 
mechanisms in GCs, and its ability to give us a unique constraint 
on the stellar type, is based on a straightforward nucleosynthetic 
prediction; if there was no Li production in the polluters, then we 
should detect positive correlations between Li and C,0, and Mg 
and hence anticorrelations between Li-N, Li-Na, Li-Al. Per con¬ 
tra, if Li production has occurred in the polluters then the common 
(anti)correlations are erased - the first- and second-generation stars 
are inclined to exhibit similar Li abund ances. This has proven to 
be the case in clusters s uch as e.g., M4 dP’Orazi & Marindl20ig; 


Mt^iarelli et alj 1201 ll ; IVillanova & Geisl^ 201T Monaco et alj 


2012h . M12 (Paper i), and M54 jMucciarelli et alj2014l) . Li produc¬ 


tion appears t o have also occur red in M5 (as we report in Paper i), 
in NG C 6397 fcind et albOO^ . and in NGC 6752 (e.g.. lShen et all 
|201(]|) . but these GCs also host a distinct extreme population that 
is somehow Li-depleted (or virtually with no Li, since only upper 
limits are available). Interestingly, the chemical behaviour of inter¬ 
mediate stars in these GCs are not the outcome of a simple mixture 
of pristine material and the composition measured in their extreme 
populations. The result implies that at least two kind of pollution 
events (thus two different kind of stellar polluters) are required to 
account for the observational constraints. In the current study we 
have identified that NGC 362 demonstrates a similar pattern to that 
found in M4 and Ml2, with no extreme population present and 
the first and second generation stars exhibiting the same Li abun¬ 
dances (within the observational uncertainties). Conversely, we de¬ 
tected a significant Li-Al anticorrelation in clusters NGC 1904 and 
NGC 2808, driven by the presence of the E population, which are 
Al-rich and Li-deficient. To account for the Li-rich intermediate 
populations in these clusters we need to postulate Li production 
between first- and intermediate second-generation stars favouring 
AGB stars as the internal polluters. We are unable to conclude if 
Li replenishment took place for the extreme cluster component and 
therefore we cannot discard any of the other polluter candidates as 
progenitors to these stars. Finally, our findings seem to suggest that 


the internal Li scatter and effective Li production is related to a 
combination of a cluster’s mass and metallicity with the relatively 
low mass and/or metal-rich clusters more efficient as far as this pro¬ 
duction is concerned. 
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